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SUMMARY 

The present paper gives a new treatment, due essen- 
tially to von Karman, of the problem of the thin airfoil. 
The standard formvlse for the angle of zero lift and zero 
moment are first developed and the analysis is then ex- 
tended to give the effect of disturbing or interference veloc- 
ities, corresponding to an arbitrary potential flow, which 
are superimposed on a normal rectilinear flow over the 
airfoil. An approximate method is presented for obtain- 
ing the velocities induced by a S-dimensional airfoil 
at a point some distance away. In certain cases this 
method has considerable advantage over the simple 
"lifting line" procedure usually adopted. The inter- 
ference effects for a S-dimensional biplane are consid- 
ered in the light of the previous analysis. The results of 
the earlier sections are then applied to the general problem 
of the interference effects for a S-dimensional biplane, 
and formulse and charts are given which permit the 
characteristics of the individual wings of an arbitrary 
biplane without sweepback or dihedral to be calculated. 
In the final section the conclusions drawn from the ap- 
plication of the theory to a considerable number of special 
cases are discussed, and curves are given illustrating cer- 
tain of these conclusions and serving as examples to in- 
dicate the nature of the agreement between the theory and 
experiment. 

L INTRODUCTION 

In the autumn of 1928 Dr. Theodor von Earman, 
in a series of lectures at the California Institute of 
Technology, presented the elements of a new approxi- 
mate theory of thin airfoils, and also gave certain ex- 
tensions and applications of the theory to the 2-di- 
mensional biplane problem. The present author was 
interested in the question of the interference effects for 
a 3-dimensional biplane and attempted the extension 
of the theory to this problem. Since the airfoil theory 
had never been published, Doctor von Karman sug- 
gested to the author that the latter work it over and 
prepare it for publication along with the biplane analy- 
sis. The following paper submitted to the National 
Advisory Committee for Aeronautics for its considera- 
tion relative to publication as a technical report, is the 
result of the effort to do this. The material in sections 



II to IV is based largely on Karman's ideas, although 
the author must accept the responsibility for the de- 
tails of the analysis, since they frequently differ widely 
from those given by Karman. In many cases, also, 
the original theory has been considerably amplified 
and extended. The 3-dimensional biplane theory itself 
has been developed entirely independently. 

H. THE THIN AIRFOIL IN AN UNDISTURBED FLOW 

The present development of the theory of thin air- 
foils, in common with many others, is based on the 
method of conform al transformation in which two 
complex planes are connected by a relation of the 
form z=f (f). Here z=x J riy is the complex variable 
for one of the planes, f=f +i? is the complex variable 
for the other plane, and/ is an analytic function of f. 
Such a relation "transforms any curve, and in partic- 
\ ular any streamline, in the f plane into a corresponding 
j curve or streamline in the z plane. If the streamlines 
j in the first plane correspond to an irrotational motion, 
the transformed flow in the other plane will have the 
same property. If the velocity of such a flow is known 
at any point in the ? plane the velocity of the trans- 
formed flow at the corresponding point in the z plane 
; is given by 

I & dz 

where g%ub* C Hpt is the appropriate component of velocity. 
Very often only the absolute magnitudes of the cor- 
responding velocities need be considered, so that the 
following simplified form of the above relation may be 
used: 



where <fr are the absolute values of the resultant 
velocities and the bars signify absolute value. For a 
very, lucid account of conformal transformation as 
applied to aerodynamics the reader is referred to 
Chapter VI of Reference 1, in which these formulae 
are deduced. 
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In the present work we shall use only the very simple 
transformation — ... 

(1) s = f+ f 

so that the transformation relation for the velocities 
has the form — 

2L 



(2) 



We start with an arbitrary flow about a unit circle 
with center at the origin of the f plane. Then (1) gives 
the corresponding flow about what may be considered 
as a straight-line airfoil in the z plane, extending 
between x=* ±2. The resultant force and moment 
acting on this straight-line airfoil can be easily deter- 
mined for any simple flow. The problem is now to 
deform this straight line into a more or less arbitrary 
airfoil shape and then determine the force and moment 
acting on this final airfoil. If y=y(x) is the equation 
of the airfoil, then by taking y as a double-valued 
function we can construct an airfoil of arbitrary camber 
and thickness. However Jeffreys has shown (Refer- 
ence 2) that for normal airfoils the effect of thickness 
is small, so that, in view of the difficulties introduced 
into the present method by the consideration of thick- 
ness, we shall confine ourselves to the discussion of 
airfoils of zero thickness. These line airfoils which 
we shall discuss may be considered as the mean camber 
lines of actual airfoils of finite thickness. We impose 
the restriction that the ordinates of the airfoils shall 
everywhere be small with respect to their chords, and 
for the purpose of actually carrying out the transfor- 
mations we also make the additional restriction that 
the leading and trailing edges shall coincide . with the 
points z=*~- 2, +2, respectively. Then y will be a 
single valued function of x and the airfoil will be a 
curve between x= ±2, which is slightly distorted from 
the original straight line. 

In accordance with (1) there will be a corresponding 
curve in the f plane which will differ slightly from the 
original unit circle. For simplicity we shall refer to 
this curve as the pseudocode. If r, 6 represent polar 
coordinates in the f plane and if on the pseudo circle 
we write r = l + e, where 6 is a variable whose value is 
everywhere small compared with 1, then the equation 
of the pseudocode may be written in complex form as 

(3) r=(i+ € y 

The equation of the airfoil is obtained from this by 
applying (1). Since e is a small quantity we neglect 
its second and higher powers and in this way obtain 
veiy simply the equation of the airfoil in complex 
form: 

z = 2 (cos d + ie sin $) 
or in the more convenient parametric form: 
3 = 2 cos 6 
y = 2e sin 0 



It will later be necessary to find the velocity along 
the surface of the airfoil, being given that along the 
pseudocode. Since both velocities will be tangential 
to the corresponding surfaces we shall require only 
the rdation between the absolute magnitudes of the 
velodties at corresponding points of the two curves, 
which is readily obtained from (2) and (3). Neglecting 
powers of « we get as a first approximation: 



(5) 



ffz 2 sin 6 



where q t and & are velodties at corresponding^ points 
P, and P f (fig. I). 
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Figure 1 

If 2fjt £# represent the indicated components of 
velodty in the £ plane and the components at 

the corresponding point in the z plane, then within 
the limits of accuracy of our approximation we may 
write the velodty at the pseudocircle as q% and that 
at the airfoil as q z . Hence, taking into account the 
conventions as to directions indicated in Figure 1, 
equation (5) becomes: 1 



(6) 



2x 2 sin 6 



(4) 



i Both equations (fi) and (0) should properly be multiplied by factor* of the form 
[1+0(«)1 where 0 («) denotes a quantity which Is of the order of magnitude of < and 
which vanishes with c Hence these two equations are exact at the unit circle and 
t are In error by quantities of the order c at tbe pseudocircle. When the relation (fl) 
! Is used, however, as for Instance In obtaining (18) from the preceding equation, It 
j Is applied either at tbe unit circle or, if at tbe pseudocircle, it Is used to transform a 
I small velocity increment which Is already of the order e. Hence the error Introduced 
into any of the expressions In which the approximate relations (fi) and (6) have been 
I employed h of the order it* in agreement with the degree of approximation throughout 
! thetheory. 
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The present method of determining airfoil char- 
acteristics is essentially the following: The conditions 
for an arbitrarily assumed flow about the straight 
line airfoil are taken as known, and the changes 
produced when the straight line is deformed into a 
curved airfoil are then calculated. The actual analysis, 
however, is carried out in the f plane and the final 
results transferred back to the z plane. Let subscripts 
zero refer to the original conditions of a flow about 
the straight line airfoil or unit circle, L e. q TQ} g# 0 
are the velocities in the f plane for the flow in which 
the unit circle is a streamline- Then in order to 
distort the flow so that the pseudocode may be a 
streamline we must superimpose on the original flow 
additional velocities q/ 9 q§'. The unique feature in 
von Karman's method is that it permits a direct 
determination of these additional velocities based on 
very simple physical reasoning. The following dis- 
cussion, while not identical with that given by von 
Karman, is the same in principle. Consider the con- 
ditions at a small element (fig. 2) in order that the 




FlGCRI 2 



pseudocircle may be a streamline, i. e. no flow across 
it. Since e is small we may, to a first approximation, 
take q% as constant along any radius between the 
circle and the pseudocircle. Then 
g r d0 = d(eg,) 

According to our notation 

2r = 2r 0 + 2r' 

2#=2»„+2»' 

but at the unit circle q ra =Q and hence 

< 7 > 

where, within the limits of our accuracy, q/ may be 
taken either at the circle or the pseudocircle, and 
similarly for q% has been omitted from the 

parenthesis of (7), since its inclusion introduces only 
terms of order £ into q/ m 



In order to find g/ we employ a known relation in 
potential theory. If <f> is a potential function, i. e. a 
scalar function whose gradient gives the velocity for 
an irrotational motion, then the general expression for 

in polar coordinates, subject to the restriction <£ = 0 
at r=c» is 

m 

v \\ A % sin nd+B n cos nB 
and the corresponding velocity components are given by 



Z>4> X~\A n sin nd+B* cos nO 

*"5r"~2j 7^ 

CD 

1 3<0 cos nd— B» sin nO 
By a little calculation it can be verified that 
2*1 h ^ cot-2-d* 

m 

SA % cos nr—B x sin nr 
-pm - 

and hence 

(8) q$ fa r) - -^J^ir ft r) cot ^ d* 

where r is an arbitrary value of the variable 0. This 
is a perfectly general result holding whenever q r and 
q$ are velocities derived from a potential function 
(i. e. for an irrotational flow) which vanishes at in- 
finity* Taking <f> as the potential function for our 
additional or superimposed flow we get at the psuedo- 
circle as a special case of (8) 

»'«■ iJ 0 V(flcot^dtf 

whicb gives, in view of (7), 

(9) ff/ (r) ± faj cot ^ 63 . 

Starting with an assumed velocity around the unit 
j circle equation (9) gives the velocity g#=ffi fl +2#' at the 
i" pseudocircle from which the velocity at the airfoil 
j may be determined from (6). In general q§ will not 
j be zero at 0=0 which implies that the velocity g a 
j at the trailing edge of the airfoil will be infinite. We 
■ now impose Kutta's condition that at the trailing 
i edge the flow must be smooth; i. e M the velocity at this 
point must be finite. Then, since 0=0 corresponds 

i It should be remarked that Tor all integrals of the form (8) the principal value 
of the Integral Is to be taken* The reader Interested In a mathematically more rigor- 
ous derivation of equation © Is referred to page Oof H. VWat's book " La Resistance 
de Fluids," Sdentia Series, GanthLer-Vfllara, Paris (1020). 
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to the trailing edge, and in view of (6), Kutta's 
condition takes the form 

2#(0)-0 

In order to satisfy this condition we must superim- 
pose upon the existing flow a circulation flow char- 
acterized by a circulation r and, by velocities at the 
airfoil or pseudocircle gr. Then Kutta's condition 
becomes 

(10) 2# (0)=»o (0)+2/ (0) + 2# r (0)-0 

This equation permits the determination of r in the 
following manner. Since y is an even function of 0 
and hence e is an odd function (cf. equation 4), there- 
fore the perimeter of the pseudocircle is to a first 
approximation the same as that of the unit circle, so 
that the tangential velocity due to a circulation r 
about the pseudocircle is, within our accuracy, the 
same as that which would exist at the unit circle due to 
a circulation r about it. Hence at the pseudocircle 

r 

2ff r (0) « constant =~— ^ 

where r is taken as positive when the tangential flow 
is clockwise about the pseudocircle or airfoil (cf . fig. 1). 
Then (1Q) gives for the determination of T. 

r-2rta, 0 (0) + 2i '(0)] 
Introducing expression (9) for g/ 

(11) r-2«a, 0 (o)- cot | de 

and the total tangential velocity at any point, 6, of 
the pseudocircle is 

(12') Ji(*)«2io(*)+a/J0)-£'<>r 

(12) a t (e)- qtt {0)-a l9 (o)-±-$*£ cot ^ dr 

+ £J?a7 ( ^ )cot i dT 

(Notice that the variable of integration has been 
changed from 0 to r.) These are the fundamental 
equations of the present airfoil theory. 

We shall first apply these equations to the deter- 
mination of the lift and moment coefficients of an airfoil 
at an angle of attack a in a uniform, rectilinear flow. 
The velocity at infinity in the z plane is taken as hav- 
ing the constant value U inclined at an angle a to the 
positive x axis. Then since the transformation (1) 
leaves the region at infinity unaltered we shall have 
these same conditions in the f plane. The velocity 



about a circle in such a flow is well known. In our 
notation we have 

j2i 0 =-2C/sin 
[2# 0 (0)=2D r sina 

In finding the lift we make use of the fortunate fact 
that a circulation is invariant to a conforms! trans- 
formation. In other words, if we are given the cir- 
culation about any closed curve in a plane, and are 
also given a second plane connected with the first by 
a conformal transformation, then the circulation about 
the corresponding curve in the second plane is identical 
with that about the original curve in the first plane. 
Hence in our case the lift of the airfoil is given by the 
familiar Kutta-Joukowsky relation 

L= P ur 

where p is the density of the fluid and r is the circulation 
about the pseudocircle as given by (11). 
Introducing expressions (13) into (11) we get 

r=47r?7sin a+2U J* l* sin <*~ a >] cofc { dd 

In simplifying this expression it is convenient to per- 
form a partial integration, but since the integral is 
improper this step requires a little investigation. Near 
8 = 0 we have along the airfoil or pseudocircle from (4) 

V «sin0 i—cxoi 

2-x 1-cos 0 t 6 0 



But near 0 = 0 ~—=*K (say) is the slope of the tan- 

gent to the airfoil at the trailing edge so that for any 
normal airfoil K is at most of order of magnitude 1. 
Hence 

- K 

for 0->O 0, .where \K\<\ 

Using this fact the ordinary methods of elementary 
calculus show that the partial integration may be 
performed and that the integrated term vanishes, 
giving 

Jo 1 — cos 6 

For all practical purposes the angle of attack is 
small, so that in the future we shall throughout call 
a a small quantity, writing 

sin a=a } cos a = l 

Hence , 

r-4T^a + 2^f^~-^(sin 0-a cos 0)d0, 
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or introducing the airfoil ordinate by means of (4) 

T=4rUa + U fVSr,- Ua Tf^U 
Jg 1 — COS V j 0 1 — cos 0 

Since y is an even function of 6 the second Integrand is 
odd and the first is even. Hence the second integral 
vanishes and we have: 



T=4tUcc+2uJ* 



cos 8 _ 

Defining the lift coefficient in the usual way, employ- 
ing however the German notation in which ^= chord, 

In our case f =4 so that 



and 

(14) 



ydd 



cos 8 



(chord =4) 



For purposes of theoretical analysis this is a very 
convenient form of the expression for Ct but for finding 
the characteristics of an. actual airfoil the following 
alternative forms are more suitable. They are obtained 
by replacing 8 by x in accordance with (4) 



(15) 



Note that in (14), (15) the airfoil chord lies along the 
x axis. 

Before discussing these equations it will be advan- 
tageous to deduce the corresponding expressions for 
the moment coefficient. Throughout the present paper 
the moment will be measured about the center of the 
airfoil — i. e., the origin in the z plane — and will be 
considered as positive when it tends to raise the lead- 
ing edge of the airfoil, i. e., stalling moments are posi- 
tive. There is, unfortunately, no simple analogue of 
the Kutta-Joukowsky equation for moments, so that 
pressures must actually be integrated over the airfoil, 
which makes this calculation somewhat more tedious 
than the corresponding one for lift. * 

The general expression for pitching moment may- 
be written 



(16) 



Al=fpx dx 



where p is the pressure at any point in the fluid and the 
path of integration is taken in a clockwise direction 



about the surface of the airfoil. Introducing the vari- 
able 8 through the substitution used before, x =2 cos 8, 

(17) Af-4 J% sin 8 cos 6 &8 

At the airfoil Bernoulli's equation gives 

where H is the total pressure head and is constant 
throughout the fluid. Substituting this in (17) the 
term containing H vanishes and replacing g z by g> 
according to (6) we have 



M 



~!J> 



p cot 8 dS 



It is convenient to consider the moment in two parts: 
Alt = moment acting on tie straight line airfoil at 
angle of attack a; M 3 = additional moment due to the 
deformation of the straight line into a # curved airfoil. 
Considering first M u *=0 for the straight line airfoil 
and hence from (12) and (13) 

2,= — 2E7[sin 0 + a(l— cos 8)] 
;.Mi= -2 pipj^lsm 2 8+a*(l -cos 8) 1 
+2 a sin 8(1 -cos 8)] cot 8 dS 



and finally 



P2r 

= +4pZ7 s o| ( 



cos^.dtf 



Mi=4xap£7 1 



As the straight line is deformed into the curved 
airfoil the velocities q* are changed by small amounts 
corresponding to the changes in g § as the circle is 
deformed into the pseudocode. We must calculate 
the additional moment Alt due to these changes. 
Let p*, g,*, q$* represent pressure and velocities for 
the straight line airfoil (y =3 €=0), and let Ap,-Ag Xr 
Lq§ represent the additional pressure and velocities 
introduced by the airfoil camber. Then from 
Bernoulli's equation 

P*+|(ffx*) J =2?*+A P -h|fe*+A 2x )* 
and neglecting the small term containing (Ag x ) 2 

Replacing p in (17) by Ap we have 

Aft = — 4p J £r*Agr r sin 8 cos 8 d8 



or 
(18) 



Mt 



9 cot 8 d8 
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g** is given by the first two terms of (12) and Aq$ by 
the last two. Hence using (13) 

Ma=4ptf 2 J%m 0+a(l~cos 9)} cotflj^J]*" ^(sinr 

-a cos t)] cot^^dr-— xt^™ T ~ a cosr)]X 
COt Jdrjdfi. 

The second integral inside the curly brackets is a 
constant with respect to 0 so that upon integration 
with respect to 6 the term containing it vanishes. 
Therefore 

M^ 2 -^ J*[co&8 + a (1-cos 9) cot ( sinT 
-a cos r)] cot -g— dTjdff—- 4 ^— J ^[€(sinr 
-a cos T )]|J^ C <5S 9 cot dfljdr 
+ a JT^ 1 * ( BinT "« cos r)][J^ (1-cos 5) X 

T — Q 1 

cot 0cot-s-~ da dr 



The second integral in the last expression for M t is 
apparently of order ae but a somewhat lengthy calcu- 
lation 6hows it to be actually of order a% so that to 
our degree of approximation it can certainly be neg- 
lected. In order to perform the integration with 
respect to 9 in the remaining part, of Ma we write 
r-9=u. Then 



x 



t — 6 

cos 0 cot — =— d0* 

0 * 



■X <■ 



COS T COS Vt 

1+COS u 



sin t sin w) — * an 

* sin ii 



so that 



T2t 

= cos t I cot ^ (1 + cos u) du 
+ sin tJ (1 + cos u) du=2ir sin t 

Ma = 4pt7* I sin r ^ [e (sin r-a cos r)] dr 
Integrating by parts 

and since e is an odd function of r 

M2=—8pTPj^ e sin r cos r dr 

or changing the variable of integration to 9 and replac- 
ing c by y from (4) 



— wX 



Combining the results for Mi and ilii we have, for the 
total moment acting on a curved airfoil at an angle of 
attack a, 

M=Mi+M 2 = 4:TrapU*-4pZpJ*y cos 8 dd 
Defining Cm by 

M=*Cx | E7V 
where in our case f =4 

(19) C*=|q<- I £y cos 0 d0 (fchord=4) 

In terms of the airfoil coordinates this gives 
n T 1 



(20) 



1 -1 f +a yadz 



ft, 



= 2 a J-iVl~^ 



(chord = 2) 



2/ cos 8 d8 



Note that in (19) and (20) the airfoil chord lies 
along the x axis. 

Expressions (15) and (20) are ossentially the same 
as those first given by Munk (Kcfcrence 3). Thoy 
illustrate very clearly that the lift coefficient of a thin 
airfoil may be split up into two parts, the first due to 
angle of attack with a constant center of pressure 
25 per cent of the chord back from the leading edge^ 
and the second due to camber whose magnitude and 
center of pressure are independent of the angle of 
attack but depend upon the camber. 

In spite of the largo amount of discussion which 
these jsquations have roceivod at the hands of various 
authors there are one or two points which should be 
mentioned explicitly here. The first relates to the 
agreoment which has been observed between these 
i expressions and experiment. The so-called angles 
of zero lift and zero moment as predicted by (15) and 
(20) are found to be verified quite satisfactorily, but 
the theoretical slope of lift and moment coefficient 
curves is not realized in practice, the discrepancies 
being roughly the same in both cases. Hence a very 
simple and satisfactory method of bringing both 
expressions into agreement with experiment is to 
multiply both by a constant factor which has been 
referred to as* the "efficiency factor" and which wo 
shall denote by 17 • 17 varies somewhat from wing 
section to wing section, but if experimental data on 
a particular section are lacking rj = 0. 875 may safely 
be taken as a good average value. 

The second point is largely one of notation. In the 
present method of derivation of the expressions for 
0 L and Cm it was necessary to take the airfoil chord 
as lying along the x axis and extending between the 
points x — —2 and +2. In the later analysis 
of section III it will be convenient to take the x axis 
in the direction of the velocity U in which case we 
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shall again require the airfoil trailing edge to lie upon 
the x axis at the point x=+2, but shall vary the 
angle of attack by rotating the airfoil chord about this 
point, so that the leading ddge will no longer be 
required to lie upon the x axis. In this case y is 
measured as before from the x axis which implies that 
for a given airfoil y=y{x) will vary with the angle 
of attack. The formulae deduced in this section remain 
valid even under this changed notation, although in 
applying them it must be remembered that a is to be 
placed equal to zero. We shall verify this fact for 
equations (14) and (19). 

Consider the same airfoil at the same angle of attack 
from the two points of view. Let primed quantities 
denote conditions relative to the new system and 
unprimed quantities conditions relative to the system 
previously employed (cf. fig. 3). 


















-2 


+2 



Fiausx 3 



Then using the old system C L and Cm are given by 
(14) and (19). In the new system we have pro- 
visionally (since we write x=2cos 6 in both cases) — 

y'de 



<V~±J> 



cos 6 
cos 6&d 



Since, however, a is small we .may write with sufficient 
accuracy 

y'={2-x) a+y = 2 (1 -cos 6) a+y 

so that 

(V=aJ^ cos* 0d0-| y cos Odd 

+ JTr?- 



or 



-cos 8 



Cm'=t> a— | y cos 6d0 



which are identical with expressions (14) and (19) for 
0 L and Cm. Hence the only restriction on the equa- 
tions of this section is that the trailing edge must lie 
on the x axis at x=+2 (x=+l for chord = 2), and 
that y must be small compared with the chord, a in 

41630—81 42 



general is then the angle between the velocity U 
and the positive x axis. The mathematical reason for 
the removal of the restriction on the leading edge 
position lies in the fact that for the convergence of the 
integrals of this section it is only necessary that y = 0 
for 0 = 0, and it is not necessary that y=Q for 6=x. 

The methods of this section may be applied to the 
calculation of the lift or moment of any portion of the 
airfoil, for example the lift and hinge moment of an 
airfoil flap, by choosing the appropriate origin and 
limi ts of integration in (16) and in the corresponding 
equation for lift. In carrying out this procedure it 
appears that very unpleasant integrals are sometimes 
encountered. In such a case an alternative method 
which avoids the* use of integrals is possible which will 
be briefly outlined here. If eg# 0 be expressed as a 
Fourier series then from (7) q/ may also be so expressed. 

q% may then be found in the same form 
by employing the series expressions pre- 
ceding (8), and finally q% may be ob- 
tained from an expression analogous to 
(12) in a form containing no integrals. 
Having q$ the determination of lift and 
moment follows the procedure already 
discussed. The forces and moments act- 
ing on an airfoil with flap have been ob- 
j tained in this manner with comparatively little diffi- 
j culty, and the earlier results of Glauert (Reference 4) 
entirely verified. 

m. THE EFFECT OF SMALL SUPERIMPOSED 
VELOCITIES 

Consider a second irrotational flow superimposed 
upon the flow of the previous section, the additional 
flow being such as would give velocities Eg? and 8q, 
at the airfoil if the latter were not present. Due to 
the principle of superposition, which states that the 
resultant of several irrotational flows is given by the 
vector addition of the velocities of the individual 
flows, the effect of th-s additional flow upon the 
airfoil characteristics can be determined by merely 
superimposing the additional velocities upon those 
which were found in Section II. It must be remem- 
bered, however, that Kutta's condition is to be 
satisfied after the superposition. 

The analysis follows closely that of the last section, 
but a few important changes must be noted. Since 
the notation tends to become cumbersome we shall not 
carry through the calculation for a general q§ Q as 
was done before, but shall immediately specify the 
particular g# 0 and j ro which are of practical interest. 
We adopt the second system of notation as described 
at the end of Section IE, taking the velocity U as 
parallel to the X axis. Since we have already seen 
that the radial velocity at the unit circle is zero for 
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the original flow without superimposed velocities, and 
from (13) q$ 0 = — 2Z7sin 0, we have 

gff 0 =« -2/7 sin 0 + 3g> 

2r 0 = *2r 

where from (5) and. {6) 

(21) 8ft=-fi&-2 sinfl 

$2r=+$2**2 sin I? 

From the two sets of equations preceding (7), the 
analogue of (7) becomes 

ff/-^l«(l7+»&) 2sin0]-5 2r 
or introducing y from (4) 

"4[»( , +*)]-* 

In order to carry out the subsequent analysis it 
will be necessary to make certain simplifying assump- 
tions as to the nature of the superimposed velocities. 
The first of these assumptions, which will be intro- 
duced at this point, is that in the cases considered in 
this paper the variations in 5q s over the airfoil chord 
are small and relatively unimportant in their effects. 
This assumption permits the variable 5q x to be re- 
placed by a constant which will be taken to be the 
value of &q z at the center of the airfoil, and will be 
written in the future simply as Then 

*'--«< ,+ *)&-* 

so that, corresponding to (9), 

In this equation and until (250 is written as ^ (r) 

to indicate that ^ is to be considered as a function of 

r not 0. From (25') on the argument (r) is omitted 
since 0 and t do not occur simultaneously and no 
ambiguity is possible. 

Introducing Kutta's condition to determine the cir- 
culation, 

r-2*to 0 (0)+ ft '(0)] 

■••s-*»*£o+*)jri-i* 



3<k are assumed to be everywhere finite, therefore 
from_.(21) 5g# (0) = Sff r (0) =>0. Hence the final equa- 
tion for the total tangential velocity at the pseudo- 
circlfi is, from (12'), " 

(22) 2 ,(0)=(^2Z7sin0-fgJ o 2 ^cot^dr 



where the terms in the first bracket are those for the. 
original undisturbed flow, and those in the second 
bracket are the additional terms introduced by the 
superimposed flow; &q v . 

In the present case it is convenient to integrate 
for both the moment and the lift. The equation (16) 
has already been given for the moment; the cor- 
responding one for the lift is 



so that 
(23) 



p sin 0d0 



in =4 p sin 0 cos 0d0 



WV follow the previous method of splitting up the 
moment (or lift) into two parts the second of which 
may be considered as a correction factor. In this 
case, however, we take the first part, L 0 or as being 
that for the curved airfoil in the undisturbed flow 
where &z* -&fr = 0. Then the additional part, AL or 
AM, is the additional effect introduced by Hq x and &q f . 
As before we let p*, q$* 9 etc., denote conditions cor- 
responding to Zo, Mq, and let Ap, Aq§ f etc., denote the 
additional pressures and velocities introduced by the 
superimposed flow, 8q v . and Cu 0 have already 
beec. found and are given from (14) and (19) by setting 
a = 0, since Z7 is parallel to the x axis. 



(24) 



Jo 1-co 

-si* 



cos 0 
cos 0d0 
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It remains to calculate AC L and AC*. Exactly as 
in Section II 3 p must be replaced byAp in (23), where 



(2#*A2i)«tf= ~^sin fljj^ ^ (r) cot ^ dr 



Hence 



2 Jo sin 5 



AM 



where £#* is given by the first bracket of (22) and Aqt 
by the second bracket. Only the portion of (g$* Ag>) 
which is an odd function of 9 will furnish any con- 
tribution to AL and AM, so that we may write 



(25) 



AAf= -p (g^Agri)^ • cot 0 



63. 



It can readily be shown that all of the integrals of 
(22) are even functions of 0. Remembering that 8 
and r are merely two different notations for the same 
variable, it follows from the geometiy of our conformal 
transformation that r= + r x and r = — r L refer at the 
airfoil to points on the top and bottom surface with 
the same value of a\ Since the airfoil is infinitely thin 
these points are infinitely close together, y, dq x , dq, 
are assumed to be continuous functions of position so 
that y (-r) = y (+r), 3£r(-r) =*£,(+ r), 5g,(-r) = 
8q r ( + r) y or y, tofr, 8q, are all even functions of r. 
From (21) it follows that 8q n &q§ are odd functions, 

and hence Sq§, 5q rj j^are all odd functions of r (or 8). 

The three integrals in (22) which contain 0 may now 
all be written in the form 

I (P) = J*f (r) cot ^ dr, where / (-r)=-/ (+r>. 
But 

7( ~* )= JT /(t) cot Sr dr 

and writing — r=<p (say) 



"Jo cot 2 dip 



But/(-*H-/(-r ?); 

/. I (- 0) - + J*/ W cot d^= I (+ ») 

Hence the three integrals containing 0 in (22) are even 
functions of 0, the other three integrals are also even 
functions since they are constants with respect to 8, 
and 8q§ is an odd function. 



-*JTtJ-H 

-n- 



cot ^ dr 



or using (21) and integrating the last term in each 
bracket by parts 

~JT^ r(r) (i+cos t) dr 

Substituting this in the expression for AL of (25) 

.-»rr^or-]*i 

and correspondingly 

AM -?^[J^I ( t ) sin r [J^cos 0 cot dfl] dr 
- (r) (1 + cos t) [J] 2 * cos 0 d0j dr 

We have seen in Section II that 

r 

Jo 



cos 0 cot 
;o 2 



? = 2t sin r 



and in the same way it is easy to show that 
J 9 2 'cot I ^d0=O 



Hence 



(25') 
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or, returning to the original variable of integration, 
6, integrating the first term of AM by parts, and 
defining AC L and AC M in the normal way for our 
chord of 4, 

ACW= -i ^?j*V cos 8to+lj*-jf (1 -cos 20) d9 

Remembering that y, Sq v are even functions of 8, and 
using (24) 

A<?i=2 S -ff C Lo + 2f* S -2f (1 + cos 6) dd 

AC M =2 S -2f 0^+~£ 5 U (1 ~ C0S 28) d6 

In order to obtain more useful expressions for the two 
integrals we express ^ as a Fourier series in 9, which 

takes the following form, since is an even function, 



«2, 

V 



COS 710 



We now introduce, our second simplification by assum- 
ing that in the cases of most interest ^ is a fairly 

slowly changing function : of position, so that at the 
airfoil it can be satisfactorily approximated by the 
first three terms of this series. We assume, therefor, 

(26) ^-aA^ + Ar cos B+A 2 cos 29 

where A^, A u A* are constants. Then 

We must now find general expressions for A, A u 
A 2 , and shall obtain such expressions in terms of an 
airfoil of arbitrary .chord t instead of for the special 
case of chord =» 4 which has previously been discussed. 
Since 9 = 0 should still represent the trailing edge .and 
B = tt the leading edge, we have as the generalization of 
(4) 

£ = | cos 9. 
Substituting this in (26) 

Letting the subscript zero denote 2 = 0 for the time 
being, 



Ay 



2\_dx\ 

A= T6[dMly)]o 

In order to simplify the notation wo shall, in the future, 

drop the subscript z&ro and let 5 ^> as well as 

represent the value of the indicated quantity at the 
center of the airfoil. Then we have . .. 

'» no 



(27) 



AC; 



SI! 



Each of the terms of (27) has a very simple physical 

jnificance. The terms containing are just those 

which would arise if the airfoil were in a rectilinear 
flow in which the velocity was increased from U to 
ZJ-EjSgx. Hence they are the increments due to the 
x component of velocity, and may conveniently be 

denoted by A t C L} A x C M . jf is the change in effective 

angle of attack at the center of the airfoil, and the 

terms in ^ are those which would appear if the angle 

of attack of the airfoil in a rectilinear flow were 

increased by Aa=^ They may be thought of as the 

terms due to the y component of velocity, and may 
therefore be written A y C Li A, C M > The term containing 

<r(^0 oecurs because °f fact that the super- 
imposed flow is such that the streamlinos at the airfoil 
have a curvature symmetrical about the airfoil conter. 
The streamline curvature has the same effect as would 
an additional camber of equal and opposite curvature 
on the airfoil in a rectilinear flow. Since this curvature 
is symmetrical about the airfoil center it furnishes no 
contribution to the moment. The term in the lift 
may be thought of as the "curvature term" and 
written A c Ck The other two terms arise from the 
fact that the superimposed streamlines have also an 
S-shaped or double curvature, and the effect is the same 
as if the airfoil in a rectilinear flow were given an 
S-shaped camber of equal and opposite amount. 
These terms may conveniently be described as the 
"double-curvature" terms and written A^Cl, A d Cv- 
All of the above physical explanations may very 
readily be verified. 

It has already been mentioned in Section II that the 
ordinary expressions for Cl and Cm can be brought 
into satisfactory agreement with experiment by mul- 
tiplying both the angle of attack and camber terms by 
an efficiency factor i?. Hence to bring (27) into 
agreement with reality the same procedure should be 
followed, i. e. all the terms due to effective angle of 
attack or effective camber should be multiplied by 1?. 
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In view of these remarks the following convenient 
method of summarizing the results of this section has 
been adopted : " 

AC l = A X C L + A,C t + A e C L + L^Cl 



(28) 



AC u = A I C u + A t C* +LiC u 



A t O L = 2rv^f 



A:C L 



where : 

C\, C^ 0 are the coefficients of an airfoil in an un- 
disturbed rectilinear flow with velocity U in the direc- 
tion of the positive z axis; 

C Lt C* are the coefficients of the same 
airfoil in the same flow but with the veloc- 
ities Bq Xi 8q, superimposed; 

Bq Al 8q, are the values of the additional 
superimposed velocities at the center of the 
airfoil chord, except when dq xy Sq 9 occur' in 
derivatives, in which case the values of the p f a 
derivatives are to be taken at the center 
of the airfoil; 

t is the airfoil chord; and 

17 is the efficiency factor which is most 
conveniently deter min ed from the fact that 
2xr is the slope (in radians) of the curve 
of ~Q L vs. angle of attack for the airfoil at 
infinite aspect ratio. 

IV. THE SUPERIMPOSED VELOCITIES FOR A 2-DIMEN- 
SIONAL BIPLANE 

In this section we shall not develop the complete 
theory of the 2-dimensional biplane, since the results 
are not of any considerable practical interest, but shall 
restrict ourselves to the determination of the disturb- 
ing velocities at an airfoil caused by the presence, in an 
otherwise rectilinear flow, of a second airfoil. The 
results so obtained will then be extended to the case 
of the 3-dimensional biplane in the next section. Since 
the results of this section are intended to be used for 
the biplane problem, an approximation method of 
finding the disturbing velocities is adopted, in which 
the disturbing velocities due to an airfoil are deter- 
mined at* a point whose distance from the airfoil 
center is of the order of magnitude of the airfoil chord 
or larger. 

In finding the velocities induced by an airfoil at a 
point some distance away the simplest and most naive 



| method consists in replacing the airfoil by a vortex 
} filament or lifting line fixed at a definite position along 
| the airfoil chord. This method, which has been widely 
used, gives the variation in induced velocity with 
variation in the lift coefficient of the airfoil as the angle 
of attack is changed. However, the effect of the varia- 
tion in the moment coefficient or center of pressure 
does not appear. In order to take account of this 
factor Prandtl, among others, has employed the device 
of t aking the vortex filament at the center of pressure 
of the airfoiL Then, as the angle of attack of the air- 
foil is varied, both the strength and position of the 
equivalent vortex filament change to correspond with 
the change in the lift and moment coefficients of the 
airfoil. This simple method suffers from the defect 
that, in finding the induced velocity at a point fixed 
with reference to the airfoil, the geometrical arrange- 
ment determined by the point, the vortex filament, 
and the undisturbed velocity U, changes with the 
angle of attack. 

Karman's method, which is here followed, takes into 
account variations both in Gl and also in G Mi but 

\ 6 Qy 




*7« 
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plane 








f ±s» - 


u u 





-Vortex 
Vortex pair 



Noti.— The vortex and vortex pair are actually superimposed. 

| without requiring any change in the geometrical 
J arrangement as these quantities vary. It consists 
j essentially in replacing the airfoil by a vortex filament 
and a vortex pair or doublet, both of which change 
in strength but not in position with variations in the 
airfoil C L and Ou. The position of both the vortex 
and the vortex pair is taken to be at the center of the 
airfoil chord. The following analysis, which justifies 
this picture and puts the results in simple analytical 
form, uses the methods of the complex potential 
function and of the elementary Cauchy theory of the 
complex variable, which are amply described in 
Chapters V to VII of Reference 1. 

Let w represent a complex potential function such 
that for the z and £ planes of Figures 1, 4 

dto 
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where f and z are connected by the original conformal 
transformation (1). Consider initiaUy a uniform, 
rectilinear flow about the straight line airfoil in the z 
plane given by 

Wo (z) = Uz 

This corresponds to the flow about the unit circle in 
the f plane; 

(29) u>o(t) = u(t + ^) 

As the straight line is deformed into a curved airfoil as 
in Figure 4 the velocity at a point P z is altered. We 
must determine the additional or induced velocities at 
P z which may be denoted by $£ r , &q v . In carrying out 
the analysis we require the airfoil to have its leading 
edge at x=— 2 and its trailing edge at x=+2 f but 
exactly as in Section II the final results will be valid 
even when the restrictions on the leading edge are 
removed. As before the analysis is largely carried 
out in the £ plane, i. e., we must find at the point Pf, 
corresponding to P^ the additional velocities intro- 
duced by the deformation of the unit circle into the 
pseudocircle. 

Let Wi be the potential function corresponding to 
these additional velocities. Then since the additional 
velocities must vanish at infinity we may write in 
general. 

where the ^4/s are complex coefficients of the form 

and a,j and bj are real quantities. 

With the scale we have chosen the coefficients A u 
A 2l etc., must have absolute values of the order of 1 
or less since, except for the circulation given by A*, 
conditions even fairly close to the circle are not vio- 
lently altered by the deformation of the circle into 
the pseudocircle. We are interested in conditions 
at P z where jR is of the order of the airfoil chord =4. 
Hence in the f plane r is of the order of 4 and the 

magnitude of the successive terms in is at most 

Ao log 4, j> pi • • ■ etc. 

Since Wi (f) is itself a small correction factor to be 
added to Wq (£) at Pf it is apparent that to any degree 
of practically interesting accuracy terms in v)i (f ) of 

higher degree than | may be neglected. Hence we 
write 

vh (f) = p(^oio g r+^) 



These terms have definite physical significances as 
follows: 

w= UAjkg (flo+ih) log £ 

represents a flow due to a source at the origin (a*) plus 
that due to a vortex at the origin (6 Q ). Since we must 
not introduce any net sources or sinks if the pseudo- 
circle is to be a closed curve, we must take a*— 0. 
Similarly 

represents the flow due to a vortex pair at the origin 
with axis perpendicular to U (a0 plus that due to a 
vortex pair at the origin with axis parallel to U (&i). 
The first gives an increase in size of the pseudocircle 
symmetrical about the x axis, and corresponds essen- 
tially to an increase in thickness of the airfoil sym- 
metrical about the mean camber line. As already 
mentioned, Jeffreys has shown that even near to the 
airfoil the effect of such a thickening is small, hence 
this effect may be entirely neglected for our purposes, 
and we may take Ai = ibi. Hence 

(30) t*(0-iZ7(Mogr + f?) 

where the first term corresponds to a vortex at the 
origin and the second to a vortex pair at the origin 
with axis parallel to U. 

The velocities in the f plane corresponding to this 
flow are given by 

and in the z plane, from the first equation of Section 

Since f is of order of magnitude 4 at Pj-, l/f a is of order 
of magnitude Kc and we have already neglected terms 
of this, order relative to 1. Hence we may take 

(cos 20-i sin 20) j 

or . 

■ 8q t , sin 6 L sin 26 

u= b «~ — bl ~l r ~ 

In order to express the velocities in terms of the z 
plane . coordinates we make another approximation 
which introduces errors of the same order as those^ 
already introduced; i. e., we take r=R and 0 = i? where 
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i?, # are polar coordinates of the point P z . Hence 
finally 

5q s , sin # L sin 2& 

0 — — 



(31) 



B 2 



5g r T cos # , t cos 2£ 



We must now find more useful expressions for the 
constants b Q , b u and shall use the well-known relations 
of Blasius for this purpose. With our conventions 
Blasius's equations for the lift and moment acting on 
a body in a flow defined by the potential function w 
are 



where the integrals are contour integrals about any 
closed curve c surrounding the body, and 3R means 
"real part of." Transforming to the {" plane 



dsoy 



/aw 



dz 
c df 



where c' is the closed curve in the f plane correspond- 
ing to c in the z plane. Writing 



/dwV 

"d^~ 
dr 



and by Cauchy's theorem on contour integrals 



•)ar 



X--|3ei(2«c 1 )-T/J(c I ) 

Af-| & (2ri [c 0 +cj]) = -xp3 (fo + cs) 
Where 3 signifies "imaginary part of." 



Applying these results to the problem in hand, we 
have from (29), (30) 

« GO =w 0 (T)+Wi GO = Tj(t+ib 0 log r+^-y^ 1 ) 
/dwy 

.. 0^1 + —- 



)(l+ i + .,..)^( 1+f9 
■■) 



l+V+2t&i , 
p +' 

/. c 1 = 2^^o^7 , 
fo+c, = -(5 0 I +2i& 1 ) Z7 1 

.*. L = 2r P lPb a 
Hence for an airfoil of chord t 

and substituting these expressions into (31) 



(32) 



fife ° L f „"n i <W\ 



sin 20 
cos 2i? 



64zi 




FiQuai 5 



It can easily be verified that these are Just the velocities 
i which would be induced by a vortex at the center of 
the airfoil whose strength was proportional to Cl, and 
a vortex pair at the same point with axis parallel to U 
and whose strength was proportional to C*. The 
procedure in this verification is exactly analogous to 
the above except that it is conducted entirely in the 
\ z plane. The vortex and vortex pair which serve to 
■ replace the airfoil are indicated in Figure 5. 
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By introducing the following notation (32) can be 
rewritten in a form more suitable for use in connection 
with the biplane problem. We shall hereafter use a 
subscript 2 to denote the upper wing and i to denote 
the lower wing of a biplane. It will often be necessary 
to use double subscripts, in which case the first sub- 
script determines the position at which a quantity is 
measured and the second indicates the cause giving 
rise to the quantity. 

For example, if Q is any such quantity, then Q 2 i is 
the value of the quantity at the upper wing which is 
produced by the presence of the lower wing. Unless 
otherwise specified, all such quantities are measured at 
the mid-point of the chord of the wing in question. 
It will be convenient to introduce the "aerodynamic 
stagger," 0, and the "aerodynamic gap," H f as indi- 
cated in Figure 5. 0 is defined as the angle between 
the perpendicular to U and the line joining the mid- 
points of the two chords, and is taken as positive when 
the upper wing is ahead of the lower. H is the dis- 
tance between the projections of the two mid-points 
on a line perpendicular to U. Since the x axis has 
been taken parallel to U the relation between /S and # 

is /3 = #-£ when the effect of the lower wing on the 

upper is considered, j8 = t?-y when the effect of the 

upper wing on the lower is considered, and that be- 
tween H and R is 

H=R cos 0. 
With these conventions equations (32) become 



(33) 



cos 3 18 cos 2/9 



-g-cos 3 P + ~4^~ \h) P P 



-^(jj) 1 cos 2 0 cos 2/3 

where the pre-subscript-^cc- has been introduced to in- 
dicate that the results are for a 2-dimensiona] biplane; 
i. e., one with infinite spans. It should be noticed 
that the second group of equations relating to condi- 
tions at the lower wing can be obtained from the first 
or upper wing group by interchanging the subscripts 
1 and 2 and replacing ft by fi + v and H by - H. This 
is a general result and permits us in the future to deter- 
mine all effects at the upper wing alone, obtaining the 
corresponding effects at the lower wing from the final 
upper wing expressions by using these simple changes. 



V. THE 3-DIMENSIONAL BIPLANE 

In this section the results of the previous sections 
are extended so as to give the lift and moment coeffi- 
cients of the individual wings of a biplane in terms of the 
assumed coefficients of the same wings when acting as 
monoplanes in an undisturbed flow. The parameters 
which enter in addition to the monoplane characteris- 
tics are the two spans, the two geometrical aspect 
ratips, the geometrical gap and stagger, the dccalago, 
and the geometrical angle of attack of the biplane 
cellule. The effects of sweepback and dihedral are 
not considered; but the latter at least should have little 
effect and can easily be handled by considering an 
equivalent biplane with no dihedral and with a gap 
equal to the mean gap of the actual biplane. In view 
of the complexity of the problem, it is unavoidable 
that the notation should become somewhat cumber- 
some, so that the various symbols and conventions 
employed are introduced in the body of the text as 
they become necessary, and the final notation is then 
summarized at the end of the paper. 

As far as the author is aware the present problem has 
been considered in a general manner only twice — by 
Betz (Reference 5) and by Eck (Reference 6) — al- 
though certain elements have been discussed by 
Prandtl, Glauert, Munk, and others. Betz uses the 
simple "lifting line" method as described in Section IV 
and the results of his theory are known to be seriously 
in disagreement with experiment. Eck does not give 
any results for biplanes of unequal spans, and his 
analysis appears to neglect one factor which may be 

of some importance. 3 Hence there seems to be some 
justification for a new consideration of the problem, 
especially in view of the considerable practical intorost 
which it holds. In the following induced drag is not 
discussed, since Prandtl's classical multiplane theory 
(Reference 7) gives the total induced drag of a biplane 
very satisfactorily. However, in the latter theory 
the distribution of lift between the two wings is assumed 
as known, so that from this point of view as well as 
from that of structural design a determination of the 
relative lifts of the two wings of a biplane is of consider- 
able importance. 

We reduce the 3-dimensional problem to an essen- 
tially 2-dimensional one in the normal manner, using 
the strip or wing element hypothesis in which the flow 
around each element of the wing along the span is 
assumed to be such that the relations of the 2-dimon- 
sional airfoil theory hold. We further assume that the 
lift and moment are uniformly distributed along the 
span of each wing whenever mutual interference effects 

1 In calculating the Induced velocity at one wing due to the trailing vortices of the 
other, Eck uses Pohlhausen's results! which are valid only for the case In which both 
wings are In the same transverse plane. Hence Eck neglects the effect of stagger on 
this portion of the downwash. 
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are under consideration. Pressure distribution exper- 
iments indicate that this assumption is not a bad one 
even in the case of biplanes with a fair amount of taper. 
It should introduce only minor errors when applied to 
such interference effects, although it is entirely unsuit- 
able for a treatment of the monoplane problem. For 
simplicity in the carrying out of the analysis both wings 
are considered as having rectangular plan forms, 
although the results are expressed in terms of span and 
aspect ratio so that they may be extended to other 
cases. The above assumptions imply that, whenever 
mutual interference effects are considered, the lift and 
moment coefficients for all elements of a particular 
wing are the same and are equal to the coefficients for 
the complete wing. The actual disturbing velocities 



strength and extending from wing tip to wing tip 
along the center of the chord. However, for the finite 
wing the two trailing vortices implied by the assump- 
tion of uniform lift distribution must also be considered. 
They are assumed to extend from the wing tips down- 
stream to infinity and to have their axes parallel to 
the velocity U. In order to fix ideas we shall^consider 
the disturbing velocities caused by the lower Jwing at 
the upper wing. Then the lower wing mustbe replaced 
by a horseshoe vortex of breadth 61 and^a vortex pair 
of the same length, as indicated in Figure Ji, where b 
denotes span. The velocity induced by^thisjsystem 
at a point P of the upper wing is to be f ound,^and then 
the mean value of this velocity over the^upperjwing. 
is to be obtained by integrating over 6 2 . 
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at one wing" due to the other vary along the span. 
Hence, in accordance with the above assumptions, they 
must be replaced by equivalent constant velocities 
obtained by. taking mean values across the span. 
When these "mean values have been deteimined they 
are to be substituted into equations (28) in place of 
the corresponding velocities which occur there. 

The first step in this procedure is to find thv dis- 
turbing velocities at one wing caused by the o'her. 
In Section IV it was shown that the disturbing veloci- 
ties due to an in finitely long wing could be calculated 
by replacing the wing by a vortex and vortex pair 
extending along the center of the wing chord. Hence 
a finite wing with constant C L} C u , and chord is to be 
replaced by a vortex and vortex pair, each of uniform 



If, as in Section IV, we write «j=the velocity 
induced at a point P by an infinitely long rectilinear 
vortex of strength r, then is the velocity which 
would be induced at P by the vortex r in a 
2-dimensional flow. If 2= the velocity at P due to a 
finite length of the vortex r, then the- well-known law 
of Biot-Savart may be written so as to give the follow- 
ing purely geometrical relation between q and «g: 



cos (p— cos y 



where <p and <p r are the angles between the vortex 
line and the lines joining the ends of the vortex seg- 
ment to P. The angles <p and <p' for the bound vortex 
of the lower wing are indicated in Figure 6. Choosing 
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the coordinate system indicated in the figure with the 
origin at the center of the lower wing, x axis parallel 
to 27, y axis perpendicular to U, z axis along the span, 
and letting r, 6, z be cylindrical coordinates of P, then 



-2- 



"0 
«2» 



2r = 0 
<o2r = 0 



2* = 



= 0 
= 0 



Hence the expression given above determines the 
relation between any component of q and the corre- 
sponding component of »£, since q and «,g are parallel. 
It might be expected that the same geometrical 
relation would also hold between the components of 
velocity induced by a finite, rectilinear, uniform 
vortex pair and those induced by the same vortex 
pair extended to infinity in both directions. An 
investigation of this point indicates that for a vortex 
pair the Biot-Savart geometrical relation does hold 
between q r and while the relation between 
q 9 and »g> is somewhat different. The actual equa- 
tion connecting q 9 and n q 9y which is too complicated 
to be conveniently used here, does not differ greatly 
from the Biot-Savart relation for cases of interest in 
the present problem. This is particularly true when 
mean values over the spans are taken. Also it must 
be remembered that the total effect due to a vortex 
pair is fairly small compared with that due to the 
corresponding vortex in most cases, here considered. 
Hence it appears that for the. purposes of this analysis 
it is satisfactory to assume that the Biot-Savart 
geometrical relation holds between the corresponding 
components of q and «g whether the velocities are 
induced by a vortex or by a vortex pair. 

(a) Effect of the Bound Vortex and Vortex Pair. 

In view of the preceding remarks and since we have 
assumed both C L and 0 M constant across the span, we 
can discuss together the induced velocity at P caused 
by the bound vortex and by the vortex pair. We have, 
therefor, for the induced velocity atP arising from both 
of these causes 



/'6gA cos <?-cos <p' 
p m\Ujn 2 

COS (p — cos <p f 
2 



where 



(^),%(l) a 
(^jy^ and (^jj) a 1 * 6 the appropriate ex- 



pressions of (33) as deduced for the two-dimensional 
biplane (infinite span). We must now find the mean 
values of these quantities over the upper span and 
shall use a bar to denote such mean values. Then 

2 

and Marly for (f)^ But and J$) 



are both constant with respect to z, whence 

&)n'XvX • kit (C08 *- C0S 9,) dx 

1% 



/a 



cos ^-cos <p') da. 



Defining the aerodynamic stagger 0 and gap II exactly 
as in the previous section, we have for the point 



cos <p— 



*+*•/> 



i coa <p = 



H 1 



A little calculation now gives, since cos *fl 3 "jy! +a i 
h _ 

i J4 (cos *' - cos ^ d2 = triiistV 1 + (^2=r cos 

a 

It is convenient to introduce the following parameters: 



(34) 



/ _ hi — I. 



- cos/S 



in terms of which we have 



Of 



'cos <p- cos tp?) dz= 



r —r 



Hence 



Now introducing the expressions given in (33) for 
cos* 0 cos 2 i\ r —t- t 



Since we are considering rectangular wings the two 
aspect ratios are defined as follows: 

A -A a =h 

M *2 
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Then the above equations may be written in the form I 



(35) 



G»+aO (r-r') sin 2/3 J 



(l).-ask*»-x x 

G*+/0 (r-r 7 ) cos2jgJ 

These expressions give the mean velocity at the upper 
wing induced by the bound vortex and vortex pair 
associated with the lower wing. 

(b) Effect of the Trailing Vortices. 

In accordance with our simplifying assumption 
which consists in replacing a wing by a horseshoe 
vortex and by a vortex pair, we have now to consider 
the mean velocities induced at the upper wing by the 
two vortices trailing downstream from the lower wing 
tips. These vortices afford no contribution to &q st 
so that all that must be calculated is their effect on 

^jj- Let us consider first the effect of the vortex 

extending from the right wing tip, z = bj2 i to infinity. 
The Biot-Savart Law gives for the resultant induced 
velocity at P 

1+COS <p 

where <p is the angle between the direction of the x 
axis and the line joining the wing tip to P. »g is 
given by 

"2=97 



where I"i is the circulation around the lower wing. 
Hence 



0" g 



2»=~2 



This velocity is perpendicular to the plane containing 
P and the vortex so that the component parallel to 
the y axis is 

6i 
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Hence 

2 r, . & 



2»' 



A ______ 

The other trailing vortex will give a corresponding 
term in which z is replaced by — z. Hence for the 
resultant vertical velocity induced at P by the trailing 
vortices: 



8r 



As before this expression must be integrated over 6 a 
in order to obtain its mean value. The integration 
gives directly 



(*).- 



Cr.fi 



27* + 



log 



+log- 



V 51 .'! _ L ')'- ^ /g,+1,+ ^) i 

By introducing the parameters defined above and 
performing a little reduction this leads to the simple 
expression 

(36) If A "CT 6_ l0g F+^P 

giving the mean velocity at the upper wing induced 
by the trailing vortices of the lower wing. 

(c) Complete Mutually Induced Velocities. 

Comb ining (35) and (36) we get for the mean value 
over the upper wing of the total velocity induced by 
the lower wing: 



I 



The relation between Ti and the lift coefficient is given 
by 

but 5<fc is small compared with U and its inclusion 
gives only second order terms in the final result, so 
that we may write with sufficient accuracy 

1 1 = — o 



In order to obtain the streamline curvature corre- 
sponding to these induced velocities we must find 
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where the values of the derivatives are taken at the 
center of the upper wing chord. From the geometry of 
Figure 6 it appears that 



d. 

dx; 
d 2 



cos a /3 d 
■ iT d> 



^^(cos^^sin^cos^) 

Remembering that p, //, r,- are functions of /?, the 
desired integrations can be carried out in a perfectly 
straightforward manner. The resulting expressions 
are a little lengthy, but the final results can be mate- 
rially simplified by the introduction of certain auxiliary 
functions which will now be defined. Consider the 



coefficient of Ci^ in (jj) 2l 



(r-r') sin j3-log 



r + sin j3 



V + sin0 

«[r sin 0-log (r+sin 0)]- [r' sin 0-log (r' + sin 0)] 
= 2<? (jx, J3) ~2e (/*', P) =*2E (say) 

Hence this quantity 2E } which is a function of //', 
0, can be expressed as the difference of two functions 
of only two variables, the form of the two functions 
being the same and only the arguments differing. 
Curves giving e in terms of /* (or /*') and P are given 
in Figure 12 from which E may readily be calculated. 
It is found that the results of. the differentiations can 
similarly be expressed in terms of the differences of 
pairs of functions, the functions in each pair having 
the same form but having arguments ji t jS and ft 
respectively. In this form the results may be sum- 
marized as follows: 



(37) 



+^0«+M')(r-r')sin2(3] 



-^(M+M'Kr-r') cos 20] 



d. 
dx 



dx* 



+^(m+m / > 32<?*] 



(38) 



sin* fi+l 



where 

G=g (m, 0 - flf (m', 0 <?* =$*(* J8) 0*', « 
and 

«G*i £)— sin/9-log (r + sin jS)] 
/^/S)=|[i + ^sin 2 /S-rcos>i8] 

=^[7 cos 20+r (6 cos a JS - 1)] 
g cos» p+ 3 £ cos 2/3 

ff* Cu, fi--j&\r (30 cos* /3-27 cos* fi+2) 
+ £(20 cos* 0-24 cos 2 0+5) 

-^(2 cos* 0-3 cos 2 0 + 1) J 

Graphs of e,/,/*, ff, ff* have been constructed and are 
given in Figures 8-12. From these curves E 9 F } F* t 
Q } G* can readily be obtained for any particular case. 

(d) Effects of Mutual Induction. 

The results given in (37) can now be introduced into 
equations (28) to give the changes in the coefficients of 
one wing caused by the presence of the other. As a 
generalisation of our former notation we write: 

Cm Ou v are fc k 0 coefficients of the indi- 
vidual wings of the biplane, and 

<?£io> "CW ^20 are tto coefficients of the indi- 
vidual wings when acting as monoplanes in an undis- 
turbed flow at the same geometrical angle of attack 
as that which they have in the biplane, when the latter 
is in the particular attitude under investigation. 
With this notation (28) becomes 



1 etc. 



A simplification is achieved by writing 

so that the changes in the upper wing coefficients due 
to the presence of the lower wing may finally bo 
written: 



(39) 
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A,Ci,=< 



A S C. 



where A w represents the total change due to mutual 
Induction. 

In order to obtain the corresponding changes at 
the lower wing the simple procedure mentioned at the 
end of Section IV may be followed; i. e., in (39) the 
subscripts i and * are interchanged, H is replaced by 
— H, and p by /S + x. From the definitions of (34) it 
follows that 

p'-*—p' 9 r— >r, r'-w' 

Similarly from (38) 

F-+F, F*-+-F*, G-*-G } G*-+G* 
E requires a little more care for 

e Oi, p)-*e 0*,-/S), « fl-w (m',-0) 
Hence if we define 

It is convenient to define 



then 



■ Making these substitutions the following expressions 
are an immediate consequence of (39) : 



(40) 



A„G 



A,C t 



Mi 4 
A^^AiC^ 

A„C Ll =(A I +A r +A e +A d )C il 
A«(7k ! = (A t + A, + Ad) ftn 

(e) Effect of Self-induction. 

To fix ideas consider first the upper wing. If we 
write 

CV.-C^+ ACL 0 Ut = Cifi, + ACv 



in which case 



then we are trying to find ACl and AC*. In general, 
Ci, differs from (7^ which implies that the downwash 
at the wing caused by its own trailing vortices is differ- 
ent in the two cases. This change in downwash causes 
a change in effective angle of attack which may be 
written as 5a a , and a corresponding change in lift and 
moment coefficients: 

£^0^=2x7150!, A t C f K,=|'j3ai=^ A.Ct, 

The subscript , indicates that the changes arise because 
of self-induction rather than from any mutual inter- 
ference effect of the two wings. The total change is 
the sum of the changes due to mutual and self-induc- 
tion, i. e. 

A<7 £a = (A,+A.) £7n, ACV,= (A« + A t ) G Ui 

In order to calculate the value of la% we must introduce 
a new assumption as to the distribution of lift along 
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the span, since the uniform distribution employed in 
calculating mutual effects leads in this case to infinite 
velocities which are physically inadmissable. Hence, 
for the present purpose we shall assume an elliptical 
lift distribution giving a constant value of da 2 along the 
span. This lack of consistency in the hypotheses 
underlying the present theory is certainly to be 
deplored from the standpoint of elegance, but it is 
probably of little importance to the actual results 
Assuming, then, an elliptical distribution 



so that 



tA 2 



or finally 
and similarly 

(410 A ^i=n^^ A A 



u 




Figuri 7 

(f) Change from Aerodynamic to Geometric Stagger 
and Gap. 

The analysis so far has been entirely in terms of the 
aerodynamic stagger 0 and gap ff, which, for a given 
cellule, vary with the geometrical angle of attack oc. 
We now introduce the " geometrical stagger" <r and 
the "geometrical gap" 6 which are defined exactly as 
in Technical ReportTNumber 240 of the N. A. C. A. 
(Reference 8), except that we take the reference point 
of each wing at the center of its chord instead of at the 
leading edge. From Figure 7 it is apparent that 



R- 



H 



G 



cos 0 cos <r 



Hence the parameters defined in (34) may be rewritten as 



bi + bi , b\ 

— - cos <r, m = 



20 



cos <r 



so that it now appears that they are constant for a 
given wing cellule and do not vary with the angle of 
attack. Hence all of the quantities in parentheses 
in (39)~and (40) are geometrical constants whoso values 
for a given biplane are independent of the angle of 
attack. If we define a as the geometrical angle of 
attack of the chord line of the upper wing (since this 
is the reference line in the determination of G and <r), 
then 

and all of the parameters involved in our problem, 
with the exception of the monoplane airfoil characteris- 
tics, are expressed in terms of geometrical constants of 
the wing cellule and the angle of attack a. 

(g) Final Results and Method of Application. 

The following form has been adopted as the most 
convenient method of summarizing the results of 
this section: 

0 L - 4- AC L &C l = L m C L + L S C L 

A*Cl =AxCt +AA +**Cl +&dC L 
Note: Cm is measured about the center of the 



(42> 



airfoil chord. 



b x + b 2 



61 — 62 



E^e Qt, 0) - e 0) E* - e Qt, - p) - e - P) 
F=j (/*, fi) -/ W, 0) F* =/*G*, 0) -/*(/, ft) 

Upper wing 



*A-^ftj)«-«»'(ti r r)«* 
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Lower wing 



^10 



sin 2/3 



i 



The method of using these results is essentially one 
of successive approximation. The following procedure 
has been found by the author to be most satisfactory: 

A series of values of a are chosen for which the 
calculations are to be made. Usually four points 
(say a=0°, 4°, 8°, 12°) will be found sufficient to 
enable continuous curves to be drawn giving the 
values of the final coefficients over the usual flying 
range. The monoplane characteristics Glq, G P. Qt 

G]£(\ 

Cu 0 , tt-' of each wing are then tabulated for the 

conditions corresponding to these angles of attack 
of the cellule. In making this tabulation it is best to 
first plot faired curves for the monoplane data and 
then pick the values used off of these curves, since 
experimental deviations from the smooth curve 
values are very much exaggerated in the succeeding 
calculations. From the monoplane data q may also 
be deter min ed or the average value of ij = 0.875 may be 
used. The quantities pf, r, r r are now calculated from 
the geometry of the wing cellule and the coefficients 

in parentheses ( ^ » v etc.^ are de- 

termined. The values of 0 corresponding to the assumed 
a's are listed and the various trigonometric functions 
required are tabulated for the a's in question. The 
magnitudes of the auxiliary functions e, — g* corre- 



sponding to these values of a (i. e. jS) are read from the 
charts and E } E* f — 6* are determined. 

Values of C^, C Mv G Ul are now assumed and 
AC Lv ACx^ kCuv &C M% are calculated. The assumed 
quantities C^, • • • Cm v may either be the monoplane 
values, 0 Lw — or they may be values estimated 
as more nearly correct by the calculator on the basis of 
his experience. In any case from the A<7 Ll , — LOu s 
so determined new values of Clv • • ■ • Gu % are obtained 
from the calculation which will in general be different 
from those originally assumed. 

The process is then repeated, using the values ob- 
tained from this first step. The results of this second 
step are then introduced and the process repeated for 
the third time. Eventually the assumed and calcu- 
lated values of G^ . . . . C M 2 will agree and these are 
then the final results. "With a little practice in as- 
suming reasonable values to begin with it is often 
possible to get the solution in a single step, and two 
steps should almost always suffice. Even when the 
first step does not lead to a solution it is often possible 
to estimate the effect of a second step and so write 
down the final result with satisfactory accuracy with- 
out actually repeating the calculations for the second 
time. 

The process is, unfortunately, somewhat lengthy and 
tedious but the author has found that, after a little 
experience, the characteristics of a biplane at four 
angles of attack can be obtained by one man in a 
comparatively few hours. A portion of the data and 
calculations for a particular biplane are given below as 
an example and should indicate fairly clearly the 
general method. 

Equal Wing Biplane 

Duta from N". A. C. A. Technical Note Number 310. [Alrfoi section— Clark Y] 

Ai=A 2 : =6 i\ = &4 

0-= +27° £ = &i/6 No decalage . 

M= 5.346 r = 5.440 i? = 0.88 0^' = . 891(7*1 
/ - 0 r' - 1.000 £=27°-a CW = .89 1 0 M i 



St 2xAi ^51 2tA 7 = 



.118 



v A 1 A 1 ~ v A i A l 



=.131 



T t ^i~hiAi~ *AiA7 ~ AzAi* -' 116 



UA± 



2yA t 



it v r-r* 5, v r-r' 

Mi 2 TtA* 2 ^°1 + 2 V /A t 1 + 2^' 

Consider a=8° 0 = 19° 



.227 



Monoplane data 



Cm 



C^ = 1.011 tf^.199 tf.P. w =30.3% ^? =.197 



CL,- 1.048 C* 10 «.199 aP. l0 =30.9% $^=.191 
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AvAoLiary junctions 
E=e (5.35, 19°) -e (0, 19°) = .02-. 03 = -.01 
E*=e (5.35,-19°)-« (0,-19°) = -1.70 -.03= -1.73 

Similarly ^=-.53 ^* = .95 <?=.26 <?* = .19 
Also cos /3 = .946 sin 2/5 = .616 cos 20=.788 

This data for a = 8° is repeated for the other angles 
of attack to be investigated. 
Equations (42) axe then rewritten in the form: 

^O^C^Cz^ (.118 cos/3) + Q Ml * C Lso (.118 sin 2/3) 
[LyG^G^ (M7E) - C Ml ' (.326 cos 20) etc. 



in which the values of all the parameters independent 
of a are introduced. Next a table of the following type 
is prepared: 



a 


.118 
cos P 


J.I8 
sin a? 


.147 
E 


.H7 
E* 


.326 
C0fl2j8 


.181 
F 


.131 
F* 


.us 


JIB 
G* 


0° 









































8°._ 


.112 


.078 


-.001 


-.-264 


.267 


-.099 


.124. 


.030 


.022 


12° 

























in which all of the values are filled in. 

Considering a particular angle of attack, «, values 
are now assumed for Cm 2 > &f V fr° m ^ e 

latter of which C'u v C'm v are calculated. Equations 
(42') are solved using these assumed values and the 
tabulated quantities. In our example (a = 8°) the fol- 
lowing is obtained: (In working this example a =12° 
was calculated first, so that in making the initial as- 
sumptions for a =8° the computer had the benefit of 
the previous results. For this reason the initial as- 
sumptions approximate more nearly to the final values 
than would be the case in general.) 
Assume: 

^ = .724 ^ = .133 ^ = .986 (7^ =.206 
Then 

A x Cl2=.091 A,C^=-.031 A,^— .065 A^ -019 

AsCAf 2 =.018 Aj,CV 3 =-.008 A d C^ 2 =-.O02 
A„(7 i2 = - .024 Aifta- .005 AC f £2 = - .019 
A m (?^ 2 = .008 A f C^ 2 =.001 AC^ 2 = .009 * 

Hence G^ = .992 Gu 2 = -208 

Similarly AC^ « .316 Atf^ = - .066 from which 

.727 , C Ml = .133 

The values of Cj^ Cu v Cu 2 originally assumed 
are now replaced by the values just determined and the 



process is repeated. In this case the final result is 
obtained after this second solution and is: 

(7^ -.992 <7 % = .725 

0*2= .208 <?*! = . 133 

After a liftle practice the calculator can often avoid 
carrying through the calculation a second time, ob- 
taining the final result by inspection of the first solu- 
tion. The present example is an instance of this. 

VI. RESULTS OP THE BIPLANE THEORY AND COM- 
PARISON WITH EXPERIMENT 

In order to investigate the nature of the results of 
the theory and to determine the agreement between 
theoretical and observed values, the characteristics 
of twenty biplanes were investigated for which ex- 
perimental data were available. These biplanes were 
all without dihedral or sweepback. In general the 
agreement between theory and experiment is reason- 
ably good, although the very considerable dispersion 
in the experimental results often makes comparison 
somewhat difficult. The sources from which experi- 
mental material was taken are given in References 
9-17, A discussion of the results for the m^re im- 
portant characteristics investigated follows. 

The results for the lift coefficient of the biplane 
cellule (C Lb ) are comparatively simple. The curves 
of Gh B when plotted against a are very nearly straight 
lines up to about # CWe, above which the theoretical 
curves begin to curve downwards. The agreement 
with experiment is in most cases perfect up to about 
% Czmax- Above this point the thooretical curves 
usually fall off somewhat more rapidly than do the 
experimental ones. Biplanes with R. A. F. 15, 
Clark Y, U. S. A. 27, and Gottingen 387 airfoil sections 
were investigated. Variations in dG L Jd<x are very 
small and are completely explained by the variations 
in the experimental values of dCiJda for the various 
monoplanes upon which the calculations were based. 
Hence we may conclude that changes in airfoil camber 
have no appreciable effect on the slope of the biplane 
lift coefficient curve. Airfoil camber does, however, 
have a small effect on the angle of zero lift of a biplane 
cellule (<xq B ). Assuming always no decalage the fol- 
lowing conclusions may be drawn: For small camber 
(R. A. F. 15) ccqb^ouc, where a^i is the angle of zero 
lift of the airfoil section acting as a monoplane. 
As the camber increases, however, c^ B tends to be- 
come lower than aw. For example, the experimental 
data used gave for the Clark Y 0^= —5.8°, while for 
an orthogonal biplane calculations (and experiment) 
gave approximately aos=— 6.4°. Similarly for a 
Gdttingen 387 orthogonal biplane the results were 
aQif=— 6.7°, o£o B = — 7.7°. Hence it appears that in- 
creasing airfoil camber tends make ot the angle of 
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zero lift of a biplane lower than that of the component 
wings acting as monoplanes. 

Staggers from —15° to +30° were investigated, 
experiment and theory agreeing in general in indi- 
cating that stagger has no appreciable influence on 
the angle of zero lift. Stagger does, however, have a 
small but noticeable effect on dC L Jda. Zero or small 
positive staggers give the lowest values, which are 
definitely below those which occur for staggers of 
+ 30°. The following table indicates the nature of 
the results for equal wing biplanes with gap-chord 
ratio of one and without decalage. The experimental 
values differ from the theoretical ones by not over 
1%, usually by much less. 



R.A.F.15 


Clark Y 


Q3tt-387 [ 




a- 




r 


dCtj/da ! 


' 0° ' 3.57 




3.30 


Q" 


3-27 | 


; 28" . 3.6C 




3.59 


30° 


,40 j 



Glauert has given a very simple method of finding 
dC LB fda for equal wing biplanes without decalage 
(Reference 18), which, is apparently quite accurate 
enough for all practical purposes, ex- 
cept in the case of negative staggers. 
In this latter case both. Glauert's and 
the present theories give values of 
C Lb somewhat higher than those deter- 
mined experimentally. In the examples 
which the author has investigated the 
present theory gives somewhat better 
agreement than does Glauert's but there 
still remains a small but apparently 
consistent discrepancy. Figure 8 gives 
curves of C Lb vs. a for a particular 
case to indicate the nature of this dis- 
crepancy. Approximately the same 
types of curves are obtained in all 
of the cases of negative stagger 
investigated. 

Decreasing the gap-chord ratio de- 
creases dC L Jda as is to be expected. 
Quantitative data on this point have 
not been obtained in this paper since the present 
theory, in view of its approximate nature, is not valid 
when the gap is much less than the chord. 

Calculations have been made for one group of equal 
wing, U. S. A. 27 biplanes at unit gap-chord ratio and 
zero stagger, but with varying amounts of decalage 
(cf. Reference 15). Changing the decalage through 
the range from — 2° to +4° had no perceptible influ- 
ence on dCi^/da, and the change in the angle of zero 
lift was exactly what would be expected if there were 
no biplane interference effects. If these results are 
valid in general they signify that Glauert's simple 
theory mentioned above may be applied to equal wing 
biplanes with as well as without decalage. 

41630—31 43 



For the effect of biplane interference on center of 
pressure position the present theory is in agreement 
with experiment in indicating that the effect is very 
small. An accurate comparison is rendered difficult 
by the very large disagreements between different sets 
of experiments. In certain cases the agreement is 
practically exact, while in others there are consider- 
able deviations of the theory from experiment. The 
following tentative general conclusions may perhaps 
be drawn. For the upper wing the theory seems to 
give quite accurate results. For the lower wing there 
appears to be a general tendency for the theoretical 
C. P/s to lie somewhat ahead of those observed, in 
certain cases as much as 4 or 5 per cent. For practical 
purposes, however, since the effect of interference is 
always small, it is probably accurate enough to assume 
that the relation between <7. P. and individual lift 
coefficient is the same in the biplane as in the case of 
a monoplane wing. 

The most interesting and important quantity dis- 
cussed in this paper is the relative lift distribution 
between the two wings. It has seemed of most prac- 
tical interest and convenience to consider this point as 
follows: Curves of the ratio "lift coefficient of the 




Fegukk 8.— K. A. F.-15 equal wing biplane 
X=B,^-L0, r— 15.ft°. 
No decalage. Experimental data from R. & M. Kos. 857 and 810. 



upper wing to lift coefficient of the biplane" have been 
plotted as ordinates against biplane lift coefficients as 

abscissae (jt^ to. Ck B ^« Unfortunately the quantity 

CzJCl b is extremely sensitive to small errors in C^ 
and Cij. For example, a very small change in effective 
decalage will change CiJO L ^ from plus to minus infin- 
ity in the neighborhood of (7^ = 0. Thus very slight 
errors in the settings of the wings, slight variations 
in profile, and very small amounts of wash-in or wash- 
out will cause large variations in the experimental 
curves, particularly in the region of small lift coeffi- 
cients. Several of the experiments were conducted by 
the pressure plotting method, in which cases it was 



660 



REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



necessary to assume the coefficient of normal force 
equal to the lift coefficient, and it was often impos- 
sible to apply any legitimate corrections for the effect 
of wind tunnel interference. In addition practically 
all of the experiments were made at comparatively 
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Figure 9.— Effect of airfoil camber on equal wing blplanea 
-4-8,^=1.0, <r=+30°. No decalage. 

small values of the Reynold's Number, so that devia- 
tions from the potential flow assumed by the theory 
might be expected to be larger than at the high Rey- 
nold's Numbers occurring in actual flight. Further- 
more , the calculations give very erratic results unless 
the monoplane experimental data are 
veiy carefully faired, since small vari- 
ations in Clj 0 and Cl^ are consid- 
erably magnified in the computation 
of interference effects. For all of the 
above reasons a very close agreement 
between theory and experiment is 
not to be expected. In certain cases 
considerable discrepancies are found, 
but on the whole the agreement is 
surprisingly good. 

The general conclusions which may 
be drawn from tie results are as 
follows: 

The effect of airfoil camber on 
Cta/^B k small but its exact nature 
has not yet been very satisfactorily 
determined. There does, however, 
appear to be a tendency for increas- 
ing camber to increase the amount 
by which OrJO LB differs from unity. 
In order to illustrate this effect, as 
well as to indicate the nature of the 
agreement between theory and exper- 
iment, curves are given in Figure 9 for three equal wing 
biplanes which are approximately identical except for 
the airfoil sections used. It should be remarked in 
connection with these curves that the Clark Y biplane 
of the series had a small effective decalage of such a 
nature that above C Lb X 0.2 the lower wing had a 



larger lift coefficient than the upper when each was 
tested separately as a monoplane. Hence for this 
biplane (above C Lb ~ 0.2) G L2 /C Lb is less than it 
would be if the two wings were actually identical and 
there were no effective decalage. For nogative 
staggers experiment and theory 
agree in substantiating the gen- 
eral conclusion stated above, 
for the cases investigated, al- 
though the effect is small. For 
zero stagger OlJO Ls differs very 
little from unity, and although 
the calculated results seem to 
confirm the above general con- 
clusion, the considerable incon- 
sistencies in the experimental 
data and the deviations of some 
of the individual results from 
those predicted by the theory 
render it difficult to make any 
generalizations in this case. 
Stagger has a very definite effect, increasing 6tagger 
in the positive direction tending to increase CzJC Lb . 
As an example of this effect, theoretical and experi- 
mental curves are plotted in Figure 10 for an R. A. F. 
15 equal wing biplane with various staggers. At 30* 
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Figure 10.— Effect of stagger on equal wing biplanes 
R. A. F.-15 biplanes .4-6,-«1.0. 

Experimental data from: K. A. C. A. T. R. 256 
«r~+30° 
«r«0° 
. andR. &M.857 

<r— 15.8° - 

stagger the British and American tests were in such 
good agreement that a single curve could be plotted to 
represent both. At zero stagger the British results 
lie considerably below the American and the theoretical 
ones and have not been included in the figure . Figures 
9 and 10 illustrate a phenomenon which has boon 
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observed in practically all of the pertinent cases 
investigated. For positive staggers the experimental 
values of C^fO^ fall off rapidly for values of Gl B} 
which are large but still well below the normal burble 
point. The theoretical curves do not exhibit such a 
behavior. The discrepancy may be a real one or it 
may be due to the low Reynold's Numbers at which 
the tests were conducted. 

The effect of decalage has been investigated with 
reference to the series of biplanes with decalage 
mentioned previously. The agreement with experi- 
ment, while not exact in every instance, is reasonably 
good, and is very much better than is the case for the 
early Betz theory. 

Decreasing the gap tends to accentuate whatever 
interference effects exist, i. e., is reduced and 
ClJCl b becomes more different from unity. In the 
neighborhood of gap-chord ratio = 1, however, the 
changes in biplane characteris- 
tics for a small change in gap- 
chord ratio are comparatively 
small. This fact furnishes a 
convenient method of deter- 
mining the range of validity of 
the present theory. For an 
R. A. F. 15 orthogonal biplane 
with gap-chord ratio =0.6 the 
theory gives values of OlJO Lb 
which are much too high, while 
for the same biplane with 
staggers of ±30° the agree- 
ment with experiment, while 
not good, is quite within reason. 
In the former case a*— 10.0, 
while in the latter it is 8.7. 
It appears from the limited 
data now available that the 
theory begins to deviate con- 
sistently from experience for 
At - greater than 7 or 8, the deviations increasing 
rapidly as ju— p' rises abov9 this value. This means 
that for normal orthogonal biplanes the theory should 
be valid for gap-chord ratios above approximately 
three-fourths. As the stagger increases either posi- 
tively or negatively the permissible gap-chord ratios 
may be somewhat reduced. This breakdown of the 
theory for small gap-chord ratios was to have been 
expected in view of the assumptions made in its 
derivation. 

As a final example the characteristics of an R. A. F. 
1 5 unequal wing biplane have been calculated. GiJClb 
and the geometrical properties of the cellule are given 
in Figure 11. It will be seen that the agreement with 
experiment is satisfactory in spite of the small gap, the 
reason being that the overhang causes /i-// to have the 
small value of 5.3. 



It should be mentioned that, following a conversa- 
tion with Prof. L. Prandtl, the theory was also devel- 
oped upon the assumption that the spacing of the 
trailing vortices is less than the span of the wingi being 
equal to the spacing between the ff rolled up" vortices 
far behind an actual rectangular wing. (Cf. Chapter 
All, Reference 1.) The formulse so obtained are not 
much more complicated than those given above, but 
the application of them to any particular, biplane 
involves considerably more labor than does the use of 
the simpler theory given here. The results appear to 
be slightly better, but it is felt that the increased 
accuracy is not enough to warrant the increase in 
complexity, so that this extension to the theory is not 
given here. 

It is perhaps interesting to note the effect of the 
relative lift distribution on the induced drag of a 
biplane cellule without decalage, as calculated from 
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Figuse II.— B. JL F.-15 unequal win* biplane 

bi ti 
» No decalage. Experimental data from R. & M. No. 097. 



PrandtTs well-known formula. The error introduced 
by assuming the lifts of the two wings proportional to 
their areas instead of taking into account the relative 
efficiencies was found to be only 3 per cent in the most 
extreme case investigated. This means that for 
"practical performance calculation it is normally 
quite permissible to neglect the effect of relative 
efficiency as is usually done by the engineer. 

vn. CONCLUSION 

The airfoil theory presented in this paper has been 
applied to the biplane problem only. It furnishes,, 
however, a general method of attacking the problem 
of the behavior of an airfoil in a disturbed flow, when-r 
ever the disturbing velocities at the airfoil are known. 
For this reason it should have a fairly extended raDge 
of usefulness. The biplane theory itself is somewhat 
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cumbersome and the application to particular examples 
involves rather tedious computations, but it is difficult 
to see at present how a much simpler theory could be 
developed for so complex a problem without the intro- 
duction of unjustifiable assumptions. The agreement 
of the present theory with experiment, as indicated 
in the examples which have been, investigated, is in 
general satisfactory, although considerable discrepan- 
cies have been found in certain cases. For the further 
investigation of these discrepancies and in view of the 
interest in the question per se, it is felt that it would be 
highly desirable to have a systematic series of experi- 
ments conducted at considerably larger Reynold's 
Numbers than any which have heretofore been 
employed in this connection. 

The author wishes to take this opportunity of ex- 
pressing his sincere appreciation of the assistance 
rendered by Mr. W. B. Oswald in checking several of 
the more lengthy differentiations and all of the numer- 
ical calculations involved in plotting the biplane 
auxiliary functions. 



Daniel Guggenheim Graduate 
School of Aeronautics, 
California Institute of Technology, 
Pasadena, Calif., March 15, 1980. 



notation employed in the final results of 
the biplane theory 

Subscript ( ) 2 implies "upper wing." 
Subscript ( )i implies "lower wing." 
C Lv Cz3= lift coefficients of the individual wings in a 
biplane 

CWy moment coefficients of the individual wings 
, . in a biplane 

measured about the center of the airfoil chord 
and defined as positive for stalling moments. 

Cxiq>~~' • C f Af 30 == characteristics of the individual wings 
when tested separately as monoplanes at the 
angle of attack corresponding to C Lv 

A x , A y , etc. =■ various components of the change in 
coefficients in passing from monoplane to bi- 
plane conditions. 

6= span. 

f = chord. 

G= geometric gap defined as in Figure 7. 
0-= geometric stagger defined as in Figure 7. 
A = aspect ratio. 

a = angle of attack of the biplane cellule defined as in 
Figure 7. 

€f — i -g* = auxiliary biplane functions given in 

" Figures 12-16. 
^efficiency factor; 27nj=slope of the monoplane lift 
coefficient curve for infinite aspect ratio. 
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Figure 14.— Biplane theory auxiliary functions 
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Figure 16.— Biplane theory auxiliary functions 
9 vs. ft, ft 
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Figubx 16.— Biplane theory auxiliary functions 



